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ABSTRACT 

Context. Our Galaxy hosts at its dynamical center Sgr A*, the closest supermassive black hole. Surprisingly, its luminosity is several 
orders of magnitude lower than the Eddington luminosity. However, the recent observations of occasional rapid X-ray flares from 
Sgr A* provide constraints on the accretion and radiation mechanisms at work close to its event horizon. 
Aims. Our aim is to investigate the flaring activity of Sgr A* and to constrain the physical properties of the X-ray flares. 
Methods. In Spring 2007, we observed Sgr A* with XMM-Newton with a total exposure of ~230ks. We have performed timing 
and spectral analysis of the new X-ray flares detected during this campaign. To study the range of flare spectral properties, in a 
consistent manner, we have also reprocessed, using the same analysis procedure and the latest calibration, archived XMM-Newton 
data of previously reported rapid flares. The dust scattering was taken into account during the spectral fitting. We also used Chandra 
archived observations of the quiescent state of Sgr A* for comparison. 

Results. On April 4, 2007, we observed for the first time within a time interval of roughly half a day, an enhanced incidence rate of X- 
ray flaring, with a bright flare followed by three flares of more moderate amplitude. The former event represents the second brightest 
X-ray flare from Sgr A* on record. This new bright flare exhibits similar light-curve shape (nearly symmetrical), duration (~3 ks) and 
spectral characteristics to the very bright flare observed in October 3, 2002 by XMM-Newton. The measured spectral parameters of 
the new bright flare, assuming an absorbed power law model taken into account dust scattering effect, are iV H =12.3*?g x 10 22 cirr 2 
and r ~2.3±0.3 calculated at the 90% confidence level. The spectral parameter fits of the sum of the three following moderate flares, 
while lower (N H =^.^*fj x 10 22 cirr 2 and T ~ 1.7^1), are compatible within the error bars with those of the bright flares. The column 
density found, for a power-law model taking into account the dust scattering, during the flares is at least two times higher than the 
value expected from the (dust) visual extinction toward Sgr A* (Ay ~ 25 mag), i.e., 4.5xl0 22 cirr 2 . However, our fitting of the Sgr A* 
quiescent spectra obtained with Chandra, for a power-law model taking into account the dust scattering, shows that an excess of 
column density is already present during the non-flaring phase. 

Conclusions. The two brightest X-ray flares observed so far from Sgr A* exhibited similar soft spectra. 
Key words. Galaxy: center, X-rays: individuals: Sgr A*, X-rays: general, radiation mechanisms: general 



1. Introduction 

Located at the center of our Galaxy, Sgr A* is the closest 
supermassiv e black hole to the solar system at a d istance of 
about 8 kpc dReidlll993h lEisenhauer etatll2003L |2005|) . Its mass 
of about 3-4xlO 6 M has been determined thanks to the mea- 
surements o f star motions (e.g. JSchodel et a l. 2002; Ghez et all 
I2003L 120051) . Amazingly, this source is much fainter than ex- 
pected from accretion onto a supermassive blac k hole. Its bolo- 
metric luminosity is on ly about 3xl0~ 9 L^dA ( Melia & Falckd 
1200 it IZhaoet alj|2003b . In particular, its 2-10keV X-ray lu- 
minos ity is only about 2 .4x1 33 erg s" 1 within a radius of 
1.5" dBaganoffetalJ 120031) . Thus, Sgr A* radiates in X-rays 
at about 11 orders of magnitude less than its corresponding 
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Eddington luminosity. This has motivated the development of 
various radiatively inefficient accretion models to explain the 
dimness of the Galactic Center black hole, e.g., A dvection- 
Dominated Ac cretion Flows (e.gjNaravan et al] 19981) . jet-disk 
models (e.g. jFalcke & M arkoff 2000), Bondi-Hoyle with inner 
Keplerian flows (e.g jMelia et al.l20 00). The recent discovery of 
X-ray flares from Sgr A* has provided new exciting perspectives 
for the understanding of the processes at work in the Galactic 
nucleus. 

The first detection of such events was found with Chandra 
in October 2000. This flare had a duration of about 10 ks, 
with a flare peak luminosity of about 1.0 x 10 35 ergs~', in 
the 2-10k eV energy range, i.e ., about 45 times the quies- 
cent state (Baganoff et alT l2"()0ll) . Further on several other X- 
ray flares were detected by XMM-Newton dGoldwurm et al.l 
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| 2003j iPorquet et alj H00l iBelanger etail 120051) and Chandra 
(Baeanoff 2003; Eckartet al. 2004, 2006, 2008; Marrone et al. 



J] 

2008). The majority of X-ray flares detected up to now have 
moderate flux amplitude with factor of about 10^45 compared to 
the quiescent state. Only one very bright flare wit h a flux ampli- 
tude o f about 160 was observed in October 2002 (Porque ret alj 
120031) . Remarkably, its peak luminosity of ~3.6xl0 35 ergs _1 
was comparable to the bolometric luminosity of Sgr A* during 
its quiescent state. Th e light curve of the X-ray flares can ex- 
hibit short (e.g., 600 s, iBaganoff et al.ll200H 200s. lPorquet et al l 
120031) but deep drops close to the flare maximum. This short- 
time scale could indicate that the X-ray emission is emitted from 
a region as small as 7/?s (~ 13J?q). 

In contrast to the near-IR (NIR) flares which appear to be 
present for up to 40% of the time (e.g., iGenzel etail 2003; 
iGhezet al.ll2004t lYusef-Zadeh etail 120061 Yusef-Zadeh et al. 
2008 in prep.) , the X-ray flaring has a much lower duty cycle, of 
typica l 1-5% dBaganoffll2003t iBelanger et alj|2005t lEckart et al l 
2006). This means either that the majority of NIR flares have no 
X-ray counterpart or that the X-ray-to-NIR ratio of some of the 
flares is too small to allow an X-ray detection above the strong, 
diffuse X-ray emission from the central parsec. When both NIR 
and X-ray flares are detected simultaneously, they show similar 
morphology in their light curves as well as no apparen t delay 
between the peaks of flare em ission (EckartetaLlE004, 2006, 
2008; Yus ef-Zadeh et al.l2006l) . The current interpretation is that 
both flares come from the same region. 

We report here the results of our Sgr A* observation cam- 
paign performed with XMM-Newton from March 30 to April 
4, 2007. The whole results of this multi-wavelength campaign 
(VLA, CSO, VLT/NACO-VISIR, HST/NICMOS, Integral) will 
be reported elsewhere (Yusef-Zadeh et al. 2008 in prep.; Dodds- 
Eden et al. 2008, in prep.). We also observed during the three 
XM M-Newton observatio ns, two bright transient sources in out- 
burst (IPorquet et alj|2007l) . The source located at about 90"-SW 
from Sgr A* has been associated with the eclipsing X-ray burster 
AX J 1745. 6-2901. Seven deep eclipses were observed as well 
as type-I bursts. The second source is located at about lO'-NW 
from Sgr A* and has been associated with the neutron star low- 
mass X-ray binary GRS 1741.9-2853 (a.k.a. AX J1745.0-2855). 
The data analysis of these two sources will be reported in forth- 
coming papers. 

We report, here, for the first time a high level of flaring ac- 
tivity with four X-ray flares, one bright and three moderate, de- 
tected in half a day. The bright flare is the second brightest X-ray 
flare detected so far from Sgr A*. 

In $2] we describe the observations and data reduction pro- 
cedure used in this work. In |3] we report the timing analysis 
of Sgr A*, and the spectral analysis of the bright flare and the 
sum of the three following moderate flares observed during this 
Spring 2007 campaign. In $H we perform a homogeneous and 
self-consistent comparison of spectral properties of the X-ray 
flares of this campaign with the reprocessed data of X-ray flares 
previously observed with XMM-Newton. Finally, in $5]we sum- 
marize our main results and discuss their possible implications. 

2. XMM-Newton observations and data reduction 

We observed Sgr A* three times with XMM-Newton in Spring 
2007 for a total exposure of -230 ks. The journal of the XMM- 
Newt on observations is g iven in Table Q] T he EPIC cam eras, 2 
MOS dTurner etalj|200lb and one pn (Stru deret aT1l2001l) . were 
operated in the full frame window mode with the medium fil- 
ter. We use the version 7.1 of the Science Analysis Software 



Table 1. XMM-Newton observation log for the Spring 2007 
campaign. 



Orbit 


ObsID 


Start time 


End time 


duration 






(UT) 


(UT) 


(ks) 


1338 


402430701 


Mai' 30, 21:29:18.1 


Mar 31, 06:28:20.8 


32.3 


1339 


402430301 


Apr 1, 15:09:03.5 


Apr 2, 19:54:45.8 


103.5 


1340 


402430401 


Apr 3, 16:40:21.5 


Apr 4, 19:47:39.2 


97.6 



(SAS) package for the data reduction and analysis, with the lat- 
est release of the Current Calibration files (CCF). The MOS 
and pn event lists were produced using the SAS tasks emchain 
and epchain. The detector light curves in the 7-15 keV energy 
range computed by these tasks show that the levels of back- 
ground proton flares were high only during the last seven hours 
of the second and the third observations, where the count rate 
exceeded the detector telemetry limit, triggering the counting 
mode. 

For the timing analysis, the contribution of the background 
proton flares was estimated using a ~3' x 3' area with a low level 
of X-ray extended-emission, located at ~ 4'-North of Sgr A* on 
the same CCD, where the X-ray emission of point sources were 
subtracted. 

To optimize in our three observations the selection of the 
events from Sgr A*, we used a two-step method based on the 
flaring activity of Sgr A* and the presence of the two bright 
transient sources. We first computed for each observation a sky 
image in the 2-10keV energy range that was used to define a 
circular region of 10"-radius centered on the X-ray source as- 
sociated with Sgr A*. From this extraction region, we built a 
preliminary background-subtracted light curvq] to identify any 
flare from Sgr A*. We identified a bright flare from this extrac- 
tion region in the last observation. Its position in each detec- 
tor was determined with the SAS task edetect_chain using 
only the flare time interval, while the positions of the other X- 
ray sources were determined using the whole exposure of the 
last observation. The angular distances between the flaring X- 
ray source associated with Sgr A* and the two bright transient 
sources (e.g., forpn 87.5" +0.2" and 599.9" +0.2", respectively) 
were then used to determine the relative position of Sgr A* in 
the two other observations. This method is independent to any 
variation of the detector position angle. Then, we use the X- 
ray c ounterparts of the Tycho-2 catalog's sources (H0 g et all 
2000) to obtain an absolute astrometry (see Appendix lAl for de- 
tails). The absolute position of the X-ray bright flare is anooo = 
17 h 45 m 40.0 s , <5j2ooo = -29°00'28.6" with a one-sigma posi- 
tional unce rtainty of l'.'O, i.e., only 075 from the radio position 
of Sgr A* (lYusef-Zadeh et all 1 19991) . Therefore, the position of 
the X-ray bright flare is fully consistent with the radio position 
of SgrA*. 

3. Results 

3.1. X-ray light curves of SgrA* 

For each observation and detector, we first built the 
source+background and the background light curves in the 2- 
lOkeV energy range with 1 s time bins starting (stopping) at the 
first (last) good time interval (GTI) of the corresponding CCD. 



1 We selected for MOS and pn the events with PATTERNS 12 and 
#XMMEA_SM, and PATTERNS 4 and FLAG=8, respectively. 
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Fig. 1. XMM-Newton/EPIC (pn+MOSl+MOS2) light curves of Sgr A* in the 2-10keV energy range obtained in Spring 2007. The 
light curves are corrected from soft-proton flaring background. The time interval used to bin the light curve is 350 s. The X-ray 
flares are labeled from 1 to 5. The horizontal lines below these labels indicate the flare durations. The quiescent level of Sgr A* 
corresponds to only 10% of the non-flaring level of these light curves (see text for details). 



We rebinned the light curves to 350 s to increase the S/N. Then, 
we subtracted the background light curve (scaled to the same 
source extraction area) from the source+background light curve. 
We scaled up, using GTI information, count rates and errors af- 
fected by the lost of exposure (e.g., due to the switch from sci- 
ence to counting mode). 

Finally, the background-subtracted light curves of the three 
detectors were summed to produce the EPIC light curves. Any 
detector missing value was inferred by the one observed by the 
other detectors using a (median) scaling factor between the de- 
tectors. 

The EPIC (pn+MOSl+MOS2) background subtracted light 
curves of Sgr A* in the 2-10keV energy range, with a time bin 
interval of 350 s, are shown in Fig.Q] During the first half of the 
exposure, the light curve is almost flat, with a non-flaring level 
of X-ray emissio n consistent with the o ne observe d with X MM- 
Newt on in 2000 dGoldwurm et al.l2003l) and 2002 dPorquet et al.l 
2003). The -50% higher non-flaring level observed in 2004 with 
XMM-Newton (Belanger et al. 2005) was du e to the contamina- 
tion by a close tran sient source in outburst dMuno etalj|2005t 
Porq uet et al.l2005l) . see also $4] Excluding Sgr A* and any tran- 
sient sources in outburst, the X-ray emission inside the 10"- 
radius region centered on Sgr A* comes mainly from one point 
source associated with the complex of stars IRS 13, the candi- 
date pulsar wind ne bula G359.95-0.0 4, and a diffuse compo- 
nent (IBaganoff et aTll2003l:lMuno et alJl2003tlWang et alj |2006). 
which contribute to about 90% of the non-flaring level in the 
2-10keV energy range. 

We identify any significant deviation from the non-flaring 
level as possible flare. A possible (#1) was observed on April 2, 
2007. On April 4, one bright flare (#2) was observed followed 
by three moderate flares (#3, #4, and #5). Figures|2]and[3]focus 
on the flare light curves with a bin time interval of 100 s that 
we use to characterize these flares. We use an iterative sigma- 
clipping of the light curve of each observation to compute the 
mean non-flaring level (dashed line) and its standard deviation. 
We define the flare time interval as the period where the light 
curve deviates from the mean non-flaring level at a confidence 
level of 95% (i.e., where the light curve is above its mean non- 
flaring level plus 1.64 times its standard deviation). Table|2]gives 



the characteristics of these X-ray flares. The weak flare #1 peaks 
just above 3 sigma, and it is, therefore, considered as reliable. 
This weak flare was independently confirmed by a simultaneous 
detection in IR by HST/NICMOS (Yusef-Zadeh at al. 2008, in 
preparation). 

The bright flare has a duration of about 3 ks, similar to the 
durati on that was observe d for the (brightest) flare of October 
2002 dPorquet et alj|2003b . Its light curve is almost symmetri- 
cal, but no significant deep drop (i.e., about 50% of flux de- 
crease) is observed in contrast to the flares of September 200 
dBaganoffetalJ 120031) and October 2002 dPorquet et alj|2003l) . 
Indeed, in the case of the October 2002 flare, the significant drop 
was observed in all three instrument light curves, while the ap- 
parent drop in the EPIC light curves of flare #2 is only observed 
in the pn light curve. However, for the flare #2 we cannot rule out 
a moderate drop in the X-ray light curve. We also computed the 
hardness ratio using the 2-5 keV and 5-10keV energy ranges, 
but we found no significant spectral change during the flare in- 
terval. 

The peak count rates of the moderate flares are 2-6 smaller 
than that of the bright flare. The durations of the moderate flares 
are 2-10 times shorter than of the bright one. The time gaps be- 
tween two consecutive flares starting from flare #2 are 5.3, 3.0 
and 1.8 hours. Therefore, four flares were observed in a time in- 
terval of only half a day. A similar group of three moderate flares 
were already observe d with XMM-Newton on 2004 March 31 
(Belanger et al. 2005), but no such preceding bright flare was 
observed. This is the first time that a such level of X-ray flar- 
ing activity from Sgr A*, both in amplitude and frequency, is 
reported. 

When the time coverage of the NIR observations allowed si- 
multaneous observations with XMM-Newton, the NIR counter- 
parts of these X-ray flares have been observed: flares #1, #4, and 
#5 were observed with HST/NICMOS (Yusef-Zadeh et al. 2008, 
in prep.); and flare #2 was observed with VLT/NACO (Dodds- 
Eden et al. 2008, in prep.). 
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Table 2. Characteristics of the X-ray flares observed in April 2007 (see Fig. |2]and[3j. 



Flare 


Day Start time (ll) 


End time <a) 


Duration Total <b) 


Peak <c) Det. <d) 


j peak (e) 


Ampl. <f) 


# 


hh:mm:ss s 


hh:mm:ss s 


s cts 


cts S -1 O" 


10 34 ergs-' 





1 


2 


15:05:03 


291913503 


15:15:03 


291914103 


600 


50.9 


0.160 


3.1 


3 3 +L2 




2 


4 


05:25:30 


292051530 


06:13:50 


292054430 


2900 


1443.0 


1.052 


20.6 


24.6+*j 


103! 2 « 


3 


4 


11:32:10 


292073530 


11:37:10 


292073830 


300 


60.4 


0.291 


5.7 


6 1 +22 

—1.2 


25 + _l 


4 


4 


14:37:10 


292084630 


14:58:50 


292085930 


1300 


212.1 


0.301 


5.9 


6 3 +23 

—1.3 


26 + _ 9 5 


5 


4 


16:45:30 


292092330 


16:58:50 


292093130 


800 


127.6 


0.426 


8.4 


8 9 +32 


37+} 3 



(a) Start and end times of the flare time interval defined as the period where the EPIC light curve deviates from the non-flaring level at a confidence 
of 95%. See text for details. 

(b) Total EPIC counts in the 2-10 keV energy band obtained during the flare interval after subtraction of the non-flaring level. 

(c) EPIC count rate in the 2-10 keV energy band at the flare peak after subtraction of the non-flaring level. 
<d) Detection level at the flare peak in cr. 

(c) 2-10 keV luminosity at the flare peak, assuming an absorbed power law model (taking account dust scatter effect), see parameter fits in Table[3] 
<f) Amplitude of the flare defined as the ratio of the 2-10 keV flare lu minosity at the peak and the 2-10 keV quiescent luminosity of Sgr A* 
observed with Chandra (i.e., 2.4 x 10 33 erg s -1 , Baganoff et al. 2003, and this work, see <j4j. 
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Fig. 2. XMM-Newton/EPIC light curve of the bright flare in 
the 2-10keV energy range. The light curve is corrected from 
soft-proton flaring background. The bin time interval is 100 s. 
The dashed and dotted lines indicate the mean non-flaring level 
and the 1.64cr threshold, respectively (see text for details). The 
dotted-dashed vertical lines indicate the time interval where the 
light curve deviates from its mean level at a confidence level of 
95%. 
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Fig. 3. XMM-Newton/EPIC background-subtracted light curve 
of the moderate flares in the 2-10keV energy range. Same defi- 
nitions, and time duration, as in Fig. [2] 



3.2. Spectral analysis of the X-ray flares 

We report here the spectral analysis of the four flares observed 
during the third observation. We used as extraction region for 
each instrument, a 10"-radius region centered on the position of 
Sgr A* determined during the bright flare time interval (§0. The 
determination of the time interval of the X-ray flare is crucial to 
prevent from any bias in the spectral analysis. For example, in 
case the flare time interval includes a significant fraction of non- 
flaring level, we find that both column density and power-law 
slope values are bias. This is more important for moderate and 
weak flares. We selected X-ray events with FLAG equal to zero, 
and since pile-up was negligible, even during the bright flare, 



with patterns 0-12 and 0^4 (single and double) for the MOS and 
pn, respectively. To extract the background spectrum we used 
the same region but limited to the non-flaring time interval. The 
response matrices and ancillary files were computed using the 
SAS tasks rmf g en and arfge n, respectively. We used XSPEC 
(version 12.4.0; Arnaud| [l996l) to fit the spectrum with X-ray 
emission models. We report in Table[3]the spectral analysis of the 
bright flare #2, and of the sum of the three following moderate 
flares (i.e., #3+#4+#5) to increase the statistics. Instead of using 
the x 1 statistic, which is not appropriate for the fitting of spectra 
with low counts (e.g., #3 +# 4+#5), we use a modi fied version of 
cstat statistic (Ca shll 1979b called the W statistic dWachter et all 
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Fig. 4. XMM-Newton EPIC spectra of the flare #2 observed dur- 
ing the April 2007 observation campaign. The spectrum for the 
non-flaring period with the same extraction region is used as 
background. Squares, triangles and stars indicate MOS1 (red), 
MOS2 (green), and pn (black) data, respectively. The lines show 
the best-fit model using an absorbed power-law model taking 
into account dust scattering (see Table lB.ll for the best fit param- 
eter values using y 2 statistic; pow model). The fit residuals are 
shown in the bottom panel. 



Table 3. Best fit parameters (using W statistic, see text for de- 
tails) of the EPIC flare spectra for absorbed pegged power-law 
(pow), bremsstrahlung (brems), and black-body (bb) models, 
taking into account dust scattering (assuming A v =25 mag). The 
errors are given at the 90% confidence level. 



Flare 


Model 




r/kX (b) 


^2-lOkeV 


C/d.o.f. 


April 2007 


#2 


pow 

brems 

bb 


12 3 +21 
10.8** 
6.6*3 


9 -1+0.3 
z " J -0.3 

6.9*' 
1.5!°;| 


16.1*11 
13.7*| 

Q 7+0.7 
y - '-0.7 


2560/2998 
2559/2998 
2562/2998 


#3+#4+#5 


pow 

brems 

bb 


8-8-3:2 
8.3*| 

4-3115 


1 7 +0 - 7 
1 - ' -0.6 

>7.0 

1 0+0.4 

'■•'-03 


5 +LS 
4 8 +26 
3 6 +0 - 7 


2117/2998 
2117/2998 
2117/2998 


October 3, 2002 




pow 

brems 

bb 


10.7*3 
6 4 +10 

"•^-0.9 


9 9+0.3 
Z - -0.3 

8.2*3 

i-6t: 


25 3 +3 <5 

21 9 +5 - 2 
iJ -* -3.6 

15 6 +0 - 9 


2728/2998 
2729/2998 
2741/2998 



(a) Nh is given in units of 10 22 cirT 2 . 
<b) The black-body temperature is given in keV. 
<c) Mean unabsorbed fluxes for the flare period in the 2-10 keV energy 



range in units of 10 12 erg cm 2 s 1 



119791) that is implemented in XSPEC (Arnaud, in prep, The 
W statistic is valid for unbinned background-subtracted spectra. 
Binning would erase information and bias the fitting result. We 
show in online Appendix|B]that for rather bright flares (e.g., flare 
#2) the parameter fits obtained using this W statistic are very 
similar to that found using the x 2 statistic (Table [BTTb ; but the 
error bars obtained using the W statistic are much better con- 
strained. We fit the spectra in the 1-10 keV energy range. We 
would like to notice that similar parameter fits are found when 
fitting in the 0.3-10keV and 2-10keV energy ranges. 

The column density (Nh) of the gas along the line-of-sight 
absor bing the X-ray photons was fitte d using the XSPEC model 
wabs dMorrison & M cCammon 1983), which uses the (old) so- 
lar abundances of I Anders & Ebiharal ( T982b . We notice that sig- 
nificant revisio ns of the solar abund ances have been made re- 
cently (see e.g., lAsplund et al.ll2005l for a review), leading to a 
decrease of the metal abundances, in particular of carbon and 
oxygen, which are the main contributors to the photoioniza- 
tion cross section above 0.3 and 0.6 keV. Therefore, using these 
revised solar abundances would incre ase the absolute d erived 
value of column density by about 50% (iGrosso et al.l2 007). with 
a negligible impact on the derived value of the photon index. We 
also include the effect of the dust scattering, which deviates soft 
X-ray photons fr om our line-of-sigh t, using P. Predehl's XSPEC 
scatter model (Predeh l & Schmittll995l) . The scatter model 
multiplies the intrinsic spectrum by exp(-r sca ), where T sca is 
the scattering optical dept h, which is given by T sca = 0.087 x 
(Ay/mag) x (E/lkeVy 2 dPredehl & Schmitt|[l995l) . From the 
/T-band extinction val ue of 2.8+0.2 mag tow ard Sgr A* obtained 
using VLT/S JJVFOM dEisenhauer et al.ll2005l) and the extinction 
law o f Riek e&LebofskvN 19851) . we inferred a dust visual ex- 
tinction value of Ay= 25.0+1.8mag. 



We fitted the spectra in turn with the following contin- 
uum models: a pegged power-law model (pegpwrli\0 model in 
XSPEC), a bremsstrahlung model, and a black-body model. The 
parameters were tied between MOS1, MOS2, and pn spectra. 
Table [3] gives for these models the best fit parameters and er- 
rors at the 90% confidence level. The best-fit parameters of flare 
#2 are rather well constrained for the three models. For the ab- 
sorbed power-law model we found /V H =12.3+^x 10 22 cm and 
r=2.3+0.3. As reported in Table lB.ll very similar best fit param- 
eters are obtained using the x 2 statistic (A ? H =12.8+ 2 ^x 10 22 cm 
and F=2. 3+0.4). For illustration purpose, we report in Fig. |4] 
the corresponding background-subtracted binned EPIC spectra 
and the data/model fit residuals obtained using the x 2 statis- 
tic. The presence of a narrow FeKa Gaussian emission lines 
(cr=10eV) is not statistically required (A^- 2 <1 for one addi- 
tional parameter). We find at 6.4 keV, 6.7 keV, and 7.0 keV upper 
limits on the equivalent widths of 154eV, 135 eV, and 145 eV, 
respectively. We also test for spectral variation index during the 
fla re. First, splitting th e spectra into two equal parts (as done 
in iPorquetetal] 120031) . we obtained /V H =13.0t 2 !?X 10 22 cirT 2 , 
r=2.5+°i and/V H =11.6;^!b< 10 22 cirT 2 , r=2. 1+0.5, for the ris- 



ing and decreasing phases, respectively. These values are con- 
sistent within the error bars. If we tied the values of the column 
density (i.e., if we assume that it is constant during the flare), 
we found F=2.4+q ^ and r=2. 2+0.4, for the rising and decreas- 
ing phases, respectively. Therefore, in both case, there is no sig- 
nificant (at a 90% confidence level) spectral variations during 
the flare, consistent with the lack of spectral variability inferred 



2 Draft available at 
ftp : //lheaftp . gsfc . nasa . gov/pub/kaa/ stat.paper . ps . 



3 We used the pegpwrlw (pegged power-law) model in which the 
unabsorbed flux over the energy range of the fit is used as normaliza- 
tion. This allows the photon index and the unabsorbed flux to be fit as 
independent parameters, and to derive directly the uncertainty on the 
unabsorbed flux. 
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above from the hardness ratio. The same conclusion is found 
when splitting the spectra into three parts: rising, top, and de- 
creasing phases. The flare #2 data are also very well fitted by 
a bremsstrahlung and black-body continuum model with a tem- 
perature of about 6.8-6.9 keV and 1 .5 keV respectively (Tables[3] 
and lB.ll ). The column density value is, as expected, dependent of 
the continuum shape. The similar column density values found 
for power-law and bremsstrahlung continua, are significantly 
higher to the value of 4.5xl0 22 cirT 2 inferred by (dust) visual 
extinction (Ay ~ 25 mag), while found only slightly higher in 
case of a black body like shape. 

The best fit parameters for the sum of the weak flares 
(#3+#4+#5) are much less constrained than those of flare #2, 
but are consistent with the brightest flares within the error bars 
(i.e., at 90% of confidence level). 

U sing a slightly hig her dust visual extinction of Ay-30mag 



Riek e et all Il989f) . as used for example in IPorquet et alj 



2003), we found very similar parameter fit values. As an exam- 



ple, for an absorbed power-law model, we found for the flare #2 



jV H =12.1+f^x 10 2Z cnr z and F=2.3±0.3 (using W statistic). 

The confidence regions of the photon index versus the unab- 
solved flux and the column density for the absorbed power law 
model including dust scattering (Ay=25 mag) will be discussed 
in the following section (see Fig. [5]). 



4. Comparison of the spectral parameters of the 
X-ray flares observed up-to-now with 
XMM-Newton 

Until now, the comparisons of X-ray flares from Sgr A* pub- 
lished in the literature have been based on spectral parameters 
obtained with different methods (e.g., different version of the 
SAS and CCF, different definition of flare time interval, scatter- 
ing effect taken into account or not, etc). Here, we compare the 
spectral properties of the April 2007 flares (#2 and #3+#4+#5) 
with those of the rapid X-ray flares observed previously with 
XMM-Newton that we obtained using exactly the same method. 



4.1. Reprocessing and spectral analysis of X-ray flares 
previously observed with XMM-Newton 

We have reprocessed the obse rvations of October 3, 2002 (or- 
bit 516, ObsID: 01 1 1350301: IPorquet et al]|2003l) and M arch 
31, 2004 (orbit 789, ObsID: 0202670601 : lBelanger et all 2005) 
with the SAS version 7.1 and the latest CCF (see Appendix let. 
We did not reprocess the X-ray flare observed on September 
4, 2001, because only a part of the beginning of the flare was 
observed dGoldwurm et a l. 2003). We did not consider here the 
X-ray flare observed on August 31, 2004 dBelanger et al.ll2005l) 
because it was contaminated by an X-ray transient in outburst, 
located at onl y 2.9"-South from Sgr A*, displayin g dips/eclipses 
at this epoch (Muno et al. 2005; Porquet et al. 2005), which pre- 
cludes the secure determination of the spectral parameters in this 
case. 

The reprocessed light curves of the October 2002 and March 
2004 flares are shown in the online Fig. IC.ll The time intervals 
of the flares have been defined using the same criteria used in 
§3.11 and reported in the online Table IC.ll As done above, the 
10"-radius extraction region has been optimized in each instru- 
ment by determining the position of Sgr A* during the time inter- 
val of the flare. The flare time interval foun d here for the Octobe r 
2002 flare is very similar to the one used in Porq uet et al.l d2003). 
While, the flare time interval for the March 2004 found here us- 



ing a time bin of 150 s is shorter ( 1200 s) than the time interval 
determined in Bel anger et alj (120051) using a much large time bin 
(~ 500 s), which contained a significant contribution of the non- 
flaring state. 

We fitted the corresponding EPIC background-subtracted 
unbinned spectra with an absorbed power law model including 
dust scattering (Ay =25 mag) and the W statistic as for the flares 
#2 and #3+#4+#5 (CL3- In the October 2002 observation, the 
non-flaring period before the flare is too short (~10ks) to allow 
a correct determination of the complex shape of the background 
spectrum of the non-flaring level, which is crucial when using W 
statistic to fit the data. Therefore, we have used the longer obser- 
vation of February 2002 (orbit 406, ObsID: 01 1 1350101) where 
no Sgr A* flare were observed. Indeed, the non-flaring spectrum 
observed on February 2002 has the same shape and flux than the 
non-flaring spectrum observed on October 2002 before the flare, 
but it much more well constrained. We notice that using the x 2 
statistic, using background spectra based either on the October 
2002 observation or on the February 2002 observation, no im- 
pact on the parameter fits are found. 

For the October 2002 flare the best fit parameters are re- 
ported in Table [3] The online Fig. IC.2I shows for illustration 
purpose the corresponding background-subtracted binned spec- 
tra and best-fit model for an absorbed power-law model using 
the x 1 statistic. As for the flare #2, fits with the x 2 statistic give 
very similar parameter fits (see Table IB. It . Although, the central 
value of the photon index, F=2.2 (+0.3), is sligh tly lower than 
the va lue, T=2.5 (±0.3), previously reported by Porque fet al.l 
(2003), both values are consistent when taking into account the 
error bars at the 90% confidence level. For this flare, both the 
position of Sgr A* du r ing th e flare and the flare time interval 
used in Porq uet et al.1 (120031) were very similar to those used 
here, and then, are not responsible for the small difference in 
the central value of the photon in dex found her e . Sinc e the first 
analysis of the brightest flare by Por quet et al.l d2003l) with the 
SAS version 5.4, the cross-calibration between MOS and pn has 
been substantially improved (see e.g., XMM-CAL-TN-0018). 
We checked that using the same ext ractions regio n , flare and 
background time intervals as used in IPorquet et al.l d2003l) and 
the new calibration used in this work, we retrieved T ~2.2. 

For the moderate flare observed on March 31, 
2004, we found jV H =4.5^xl0 22 cnr 2 and F=0.8+0.6 
(C/d.o.f. = 1860/2998). Similar values are obtained with 
the x 2 statistic (jV H =4.7+:^xl0 22 cnr 2 and r=0.9+J°, 
^ 2 /d.o.f. = 18.4/21; see Fig. IC.3b . This central photon index 
value is smaller than the one reported by Belan ger et alj (2005) 
(r=1.5jg with the error bars calculated at the 68% confidence 
level), though compatible within the large error bars. The 
difference is mainly due to the method used here: a better deter- 
mination of the flare time interval, which helps to decrease the 
contribution of the non-flaring level; a spectral fit done without 
binning the spectra when using the W statistic, in order to not 
loose spectral information and then to not bias the parameter 
fit result; the dust scattering effect taken into account. A closer 
look to the fit of the March 2004 spectra show that when the 
instruments are fitted separately, very different central fit values 
are obtained (although compatible within the very large error 
bars). Moreover, the inferred column density value is lower than 
the one determined during the quiescent state (see below) and 
therefore appears unphysically low. When the column density 
is fixed to the value found for the brightest flares, the spectral 
index is 1.8+0.4 using the W statistic (C/d.o.f= 1872/2999) 
and 2.3+0.6 using the x 2 statistic (y 2 /d.o.f.=24/22), which is 
consistent within the error bars with the spectral index value 
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Fig. 5. Confidence regions of the spectral parameters of the X-ray flares from Sgr A* for an absorbed power-law model taking into 
account the dust scattering. Left panel: photon index versus unabsorbed flux in the 2-10keV energy range. The dashed, continuous, 
and dotted contour levels correspond to confidence levels for two interesting parameters of 68%, 90%, 99%, respectively (i.e., to 
AC = 2.3, 4.61, 9.21, respectively). The confidence regions corresponding to the quiescent state of Sgr A* are inferred from the 
spectral analysis of four archived Chandra observations, using the same spectral models (red contour levels) or including a Gaussian 
emission line (blue contour levels). See text for details. Right panel: photon index versus column density. The dashed-dotted vertical 
line indicates the colum n density (4.5 xlO 22 cm' 2 ) derived from the /T-band ex tinction value towa r d Sgr A * (2.8 mag with one-sigma 
un certainty of 0.2 ma g; Eisenhau er et al.l 12005) using the extinction law of Rie ke & Lebofskvl (1 19851) and the empirical relation 
of iPredehl & Schmittl (Il995l) : the dashed, continuous, and dotted vertical lines correspond to 68%, 90%, 99% confidence levels, 
respectively. 



of the brightest flares. Therefore, we conclude that the spectral 
properties of this moderate flare are poorly constrained, and we 
will not considered it in our comparison of the flare properties. 

4.2. Confidence regions of the spectral parameters 

We show in the left and right panels of Fig. [5] the contour plot 
of the confidence regions of the photon index versus the unab- 
sorbed flux and the column density, respectively, at the confi- 
dence levels of 68%, 90%, and 99% (corresponding to AC=2.3, 
4.61, and 9.21, respectively, for two interesting parameters). The 
fluxes of the flares span a large range, i.e., about a factor 5. The 
flare #2 has a flux lower than the October 2002 flare at the 90% 
confidence level, but has similar best fit values of column den- 
sity and spectral index. Due to a lower S/N, the physical pa- 
rameters of the sum of the moderate flares (#3+#4+#5) are less 
constrained, but has column density and spectral index values 
compatible (within the error bars) with those of the bright flares. 

The column density values during the flares are at least two 
times higher than the value expected from the 25 mag of (dust) 
visual extinction toward Sgr A*, i.e., 4.5xl0 22 cirT 2 . To know 
whether this excess of (gas) column density is related to the 
flare phenomena itself, we need to estimate the column den- 
sity of the quiescent spectrum of Sgr A* with the same spec- 
tral model. However, the quiescent emission of Sgr A* can only 
be resolved with Chandra. The spectral properties of the quies - 
cent emission of Sgr A* was reported by Bagan off et alj (12003). 
based on the first Chandra observations of the Galactic Center 



obtained in 1999 with a total effective exposure of about 41 ks. 
More Chandra observations of Sgr A* are now available in the 
archives. Starting from the level 2 of the archive event lists, we 
used CIA0 4 . 8 and CALDB 3.4.2 to produce Chandra X-ray 
light curves of Sgr Afl We selected four Chandra observations 
where Sgr A* displayed no flare^]. We extracted the correspond- 
ing qu iescent spectra following the method of Baganoff et al. 
(120031) : we refined the background 10"-radius extraction region 
by excluding the candidate pulsar wind nebula G359.95-0.04. 

We have performed the simultaneous spectral analysis 
of these four quiescent spectra of Sgr A* obtained with 
Chandra. First, we fit the data with an absorbed power- 
law model using the x 1 statist ic and a spec tr al bin ning of 
10 counts per bin as done in Bag anoff et al.l (12003b . When 
the dust scattering is not taken into account, we found 
A r H=9.7+|-^x 10 22 cirT 2 and F=2.5^ Ct 2 /d.o.f.=143/130) and 
an unabsorbed 2-10keV flux of 2.8 + [j 10' 13 erg cm" 2 s~' 
These values are similar to those reported by Bag anoff et al.l 
d2003l) but are better constrained. An excess of emission 
is seen near 6.5 keV, therefore we add a Gaussian emis- 
sion fine and found /V H =ll.l! 2 Jxl0 22 cm- 2 and F=3.0+^ and 
an unabsorbed 2-10keV flux of 3.2+^xl0~ 13 ergcirT 2 s _1 , 
E lm e=6.61!° : !« keV, cr line =0. 15!°;j° keV, and EW line = l.l^ keV 
(y 2 /d.o.f. = 123/127). As found by iBaganoff et all <|2003), we 



4 See Science threads of CIAO 4 . Q 

5 Namely ObsID 3665, 4683, 5950, and 5951 with exposure times of 
about 90, 50, 48, and 44 ks, respectively. 
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Fig. 6. Chandra quiescent spectra of Sgr A*. Black, red, green, 
and blue data points correspond to ObsID 3665, 4683, 5950, and 
5951, respectively. The lines show the best-fit model (y 2 statis- 
tic) using an absorbed power-law model taking into account the 
dust scattering. The fit residuals are shown in the bottom panel. 

found a steeper power-law slope when an emission line is 
added to the fit. In this work, thanks to a much longer ex- 
posure time, we are able to constrain the line energy, which 
is consistent with a He-like iron line. For illustration purpose, 
we show in Fig. [6] the background-subtracted binned Chandra 
spectra and the data/model fit residuals obtained for an ab- 
sorbed power-law model taking into account the dust scatter- 
ing effect with the x 1 statistic. To perform a homogeneous 
comparison with the spectral analysis of the XMM-Newton 
X-ray flares reported in this work, we then fit the unbinned 
Chandra quiescent spectra using the W statistic and we took 
into account dust scattering effect. We found A^h=8.4^]^x 
10 22 cm" 2 , r=2.4+°-4 (C/d.o.f.=1825/2457), corresponding to an 
unabsorbed 2-10keV flux of 3.3^;jxl(r 13 ergcirr 2 s and a 
luminosity of 2.4 + " ^xlO 33 e rg s" 1 (assumin g d-8 kpc), simi- 
lar to the value reported by Bagano ff et al] (|2003). The con- 
tour plots corresponding to this fit are reported on Fig. [5] 
If we add a Gaussian emission line we found A?h=9.8!|^x 
10 22 cirr 2 , F=3.0+0.5, E toe =6.67^;^keV, cr line <0.19keV, and 
EW line =0.9^ A 3 keV (C/d.o.f.=1781/2454). As shown in Fig. 
the column density and power-law index values found for the 
quiescent state (Gaussian line added or not to the fit), are com- 
patible, within the error bars calculated at 90% confidence level 
for two interesting parameters, with those found for the bright 
flares. 

5. Summary and discussions 

Here, we reported the data analysis of the XMM-Newton April 
2007 campaign observation of Sgr A* (three observations with 
a total of exposure of ~230ks). We observed four X-ray flares 
occurring within only half a day on April 4: one bright flare (#2, 
with a detection level of ~21 cr) followed by three moderate ones 
(#3, #4, and #5, with a detection level of ~6<t, ~6<t, and ~8<x, 
respectively). As well a weak flare (#1, with a detection level of 
~3cr) was observed on April 2, 2007. The flare #2 is the second 
brightest X-ray flare detected so far from Sgr A* and has a sim- 
ilar duration (~3 ks) to that of the October 2002 brightest flare 
(Porq uet et al]|2003l) . Its light curve is almost symmetrical but 



without no significant short-time scale drop (i.e., about 50% of 
flux decrease) contra ry to that reported for the Octo ber 2000 flare 
dBaganoffet al.ll200lh and the October 2002 flare dPorquet et al] 
2003). However, for the flare #2 we cannot rule out a moderate 
drop in the X-ray light curve. A similar group of three moderate 
flar es were already obser ved with XMM-Newton on 2004 March 
31 dBelanger et al. 12 005). but no such preceding bright flare was 
observed. 

This is the first time that a such level of X-ray flaring activ- 
ity from Sgr A*, both in amplitude and frequency, is reported. 
Observations such as those reported in this paper can eventually 
constrain flare models in Sgr A*, perhaps even ruling some out. 
For example, the quick succession of several events separated by 
only a few hours, might argue against a disruption mechanism 
that relies on the temporary storage of mass and energy, if all the 
corresponding accretion energy is released during the flare. The 
accretion rate M in th is system (so me estimates place it as low as 
10 16 - 10 17 g s _1 ; see lMelia]|2007l) might not produce a transient 
accumulation of mass AM between flares of sufficient magni- 
tude for T]GMAM/3Rs to provide the observed outburst power. 
In this expression, Rs is the Schwarzschild radius, and r] is the 
(poorly known) radiative efficiency in Sgr A*, believed to be at 
most a few percent. If instead the flares are due to a magneto- 
rotational instability, then the energy liberated during the flare 
must also be accumulated over the short inter-burst period. The 
low M, from which the energy is derived, might argue against 
this type of mechanism as well. 

On the other h and, if the flare arises from the infall of a 
clump of gas (e.g., iLiu et al] 120061; iTagger & Melial l2006) then 
th ere would be l ess re striction on how often these could come 
in. iGenzel et al] (120031) have shown that the total energy release 
> 10 39 ' 5 erg during a flare requires a gas accreted mass of a 
few times 10 19 g (assuming a radiation efficiency of ~10%), i.e. 
comparable to tha t of a comet or a small asteroid. Recently, 
ICadez et al] Q2006) have argued that the flares could be pro- 
duced by tidal captures and disruptions of such small bodies. 
The comet/asteroid/planetesimal idea for depositing the addi- 
tional mass and energy to initiate a flare is attractive. The dis- 
tance from Sgr A* at which such a small body would get tidally 
disrupted is the Roche radius, which is for a rigid body: 

*u.3.2x( *f r x (-&-,r. 

R s \4xlO 6 M / \lgcm- 3 / 

where Mbh is the black hole mass and p p is the density of the 
rigid body. Thus, for a black-hole mass of 4 x 10 6 M Q and a den- 
sity of 1 gcrrT 3 , this corresponds to 13.2 Rs, in good agreement 
with the size of the region where the flares are thought to occur. 
The flaring rates would then depend on processes occurring 
much farther out, so the fact that so many flares are seen so close 
together on some days, and much less frequently at other times, 
would simply be due to stochastic events. However, it would 
still be difficult to distinguish between a compact emissio n 
region and emission within a jet (see, e.g., Markoff et "aT]|200l1) . 
since these disruption events could still end up producing an 
ejection of plasma associated with the flare itself. 

Four of the five X-ray flares observed during this campaign 
have a simultaneous NIR flare counterpart: flares #1, #4, and 
#5 were observed with HST/NICMOS (Yusef-Zadeh et al. 2008, 
in prep.); and flare #2 was observed with VLT/NACO (Dodds- 
Eden et al. 2008, in prep.). The flare #3 has not been covered 
by a simultaneous NIR observation. This strengthens the re- 
lationship observed up to now between X-ray and NIR flares 
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when there is a simultaneous X-ray/NIR observation coverage: 
all X-ray flares have an NIR flare counterpart, while all NIR 
flares are n ot each time associated with an X-ray flare coun- 
terpart (e.g., lEckart et al.ll2004 120061: lYusef-Zadeh e7all l2006t 
iHornstein et alj|2007t iMarrone et al.ll2008h . 

We have made the first detailed comparison of X-ray flare 
properties observed with XMM-Newton (this April 2007 cam- 
paign and previous observations) based on a fully self-consistent 
analysis approach, namely we use the same SAS version with the 
latest calibration, the same definition of the flare interval (pre- 
venting from significant contribution of the non-flaring level), 
an optimized extraction region for each instrument, and we took 
into account the dust scattering effect. We found that the phys- 
ical parameters of the flare #2, iVH=12.3*jgX 10 22 crrT 2 and 
r=2.3+0.3, are well constrained and are very similar to that of 
the October 2002 (brightest) flare, N H =12.3+}!?xl0 22 cirT 2 and 
T=2.2+0.3. The spectral parameter fit values of the sum of the 
three moderate flares following the bright flare observed on April 
4, 2007, A^ H =8.8+^x 10 22 cnr 2 and T=1.7+^, while lower, are 
compatible within the error bars with those of the bright flares. 
The column density found during the bright flares is at least 
two times higher than the value expected fro m the (dust) vi- 
sual e xtinction toward Sgr A* (Ay ~ 25 mag; Eisenhaue r et alj 
120051) . i.e., 4.5xl0 22 cirT 2 . However, our fitting of the Sgr A* 
quiescent spectra obtained using four Chandra observations for 
a total observation time of ~ 230 ks, taking into account the 
dust scattering, shows that a column density excess is already 
present during the non -flaring phase. One possible explanation, 
as already pinpoint by Maed a et al.l d2002l) . is that such column 
density excess in the line-of-sight of Sgr A* could be due to 
the warm (dust free) gas associated with the ionized gas halo, 
Sgr A West extended, suggested by the turnover abs orption ob- 
serve d at 90 cm, embedded th e Sgr A complex (e.g., Pedlar et all 
1989t lAnantharamaiah etail 119991: lYusef-Zadeh et all l2000h . 
Ananfharamaiah et al.l ([T999) inferred that this ionized halo has 
a dimension of about 4' (~9pc) and an electron density of about 



100-1 000 crrT 3 . For a column excess of about 2.9xl0 22 crrT 2 as 
found for the Sgr A* quiescent state, this would correspond to a 
consistent electron density of about 1300 cm -3 . However, in case 
the genuine continuum shape is curved compared to a power law 
and bremsstrahlung models, as for example a black body like 
shape, no significant column excess would be required. 

On the basis of our current analysis we conclude that the 
dichotomy of the spectral index between moderate and bright 
flares (i.e. , a correlation betwee n photon index and flux) sug- 
gested by Belanger et alJ d2005h is not confirmed. However, to 
test this more thoroughly we need further observations so as 
to extend the sampling of the flaring activity seen from Sgr A* 
(and in particular those instances where good signal to noise can 
be achieved through the co-adding of weak to moderate flaring 
events). 

In conclusion, this study establishes that the two brightest 
X-ray flares observed so far from Sgr A* exhibited similar (well 
constrained) soft spectra. Therefore, any model proposed to ex- 
plain the flaring behavior of Sgr A* must take into account this 
observational X-ray spectral property. 

Acknowledgements. The XMM-Newton project is an ESA Science Mission with 
instruments and contributions directly funded by ESA Member States and the 
USA (NASA). 



the Galaxy, ed. H. Falcke, A. Cotera, W. J. Duschl, F. Melia, & M. J. Rieke, 
422 

Anders, E. & Ebihara, M. 1982, Geochim. Cosmochim. Acta, 46, 2363 
Arnaud, K. A. 1996, in ASP Conf. Ser. 101: Astronomical Data Analysis 

Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17 
Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Astronomical Society of 

the Pacific Conference Series, Vol. 336, Cosmic Abundances as Records of 

Stellar Evolution and Nucleosynthesis, ed. T. G. Barnes, III & F. N. Bash, 25 
Baganoff, F. K. 2003, in Bulletin of the American Astronomical Society, Vol. 35, 

Bulletin of the American Astronomical Society, 606 
Baganoff, F. K., Bautz, M. W., Brandt, W. N., et al. 2001, Nature, 413, 45 
Baganoff, F. K., Maeda, Y., Morris, M., et al. 2003, ApJ, 591, 891 
Belanger, G, Goldwurm, A., Melia, F, et al. 2005, ApJ, 635, 1095 
Cadez, A., Calvani, M., Gomboc, A., & Kostic, U. 2006, in American Institute 

of Physics Conference Series, Vol. 861, Albert Einstein Century International 

Conference, 566-571, [arXiv:0708.2696] 
Cash, W. 1979, ApJ, 228, 939 

Eckart, A., Baganoff, F. K., Morris, M., et al. 2004, A&A, 427, 1 
Eckart, A., Baganoff, F. K., Schodel, R., et al. 2006, A&A, 450, 535 
Eckart, A., Baganoff, F. K., Zamaninasab, M., et al. 2008, A&A, 479, 625 
Eisenhauer, F, Genzel, R., Alexander, T., et al. 2005, ApJ, 628, 246 
Eisenhauer, F, Schodel, R., Genzel, R., et al. 2003, ApJ, 597, L121 
Falcke, H. & Markoff, S. 2000, A&A, 362, 113 
Genzel, R., Schodel, R., Ott, T., et al. 2003, Nature, 425, 934 
Ghez, A. M., Duchene, G., Matthews, K., et al. 2003, ApJ, 586, L127 
Ghez, A. M., Salim, S., Hornstein, S. D., et al. 2005, ApJ, 620, 744 
Ghez, A. M., Wright, S. A., Matthews, K., et al. 2004, ApJ, 601, L159 
Goldwurm, A., Brion, E., Goldoni, P., et al. 2003, ApJ, 584, 751 
Grosso, N., Bouvier, J., Montmerle, T., et al. 2007, A&A, 475, 607 
H0g, E., Fabricius, C, Makarov, V. V, et al. 2000, A&A, 355, L27 
Hornstein, S. D., Matthews, K., Ghez, A. M., et al. 2007, ApJ, 667, 900 
Liu, S., Melia, F, & Petrosian, V. 2006, ApJ, 636, 798 
Maeda, Y, Baganoff, F. K., Feigelson, E. D., et al. 2002, ApJ, 570, 671 
Markoff, S., Falcke, H., Yuan, F, & Biermann, P. L. 2001, A&A, 379, L13 
Marrone, D. P., Baganoff, F. K., Morris, M., et al. 2008, ApJ, in press, 
[arXiv:0712.2877] 

Melia, F. 2007, The Galactic supermassive black hole (The Galactic supermas- 

sive black hole, by Fulvio Melia. Princeton, NJ: Princeton University Press, 

2007, ISBN 0691095353.) 
Melia, F. & Falcke, H. 2001, ARA&A, 39, 309 
Melia, F, Liu, S., & Coker, R. 2000, ApJ, 545, LI 17 
Morrison, R. & McCammon, D. 1983, ApJ, 270, 119 
Muno, M. P., Baganoff, F. K., Bautz, M. W., et al. 2003, ApJ, 589, 225 
Muno, M. P., Lu, J. R., Baganoff, F. K., et al. 2005, ApJ, 633, 228 
Narayan, R., Mahadevan, R., Grindlay, J. E., Popham, R. G., & Gammie, C. 

1998, ApJ, 492, 554 
Pedlar, A., Anantharamaiah, K. R., Ekers, R. D., et al. 1989, ApJ, 342, 769 
Porquet, D., Grosso, N., Belanger, G., et al. 2005, A&A, 443, 571 
Porquet, D., Grosso, N., Goldwurm, A., et al. 2007, The Astronomer's Telegram, 

1058, 1 

Porquet, D., Predehl, P., Aschenbach, B., et al. 2003, A&A, 407, L17 

Predehl, P. & Schmitt, J. H. M. M. 1995, A&A, 293, 889 

Reid, M. J. 1993, ARA&A, 31, 345 

Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 

Rieke, G. H., Rieke, M. J., & Paul, A. E. 1989, ApJ, 336, 752 

Schodel, R., Ott, T., Genzel, R., et al. 2002, Nature, 419, 694 

Striider, L., Brief U, Dennerl, K., et al. 2001, A&A, 365, L18 

Tagger, M. & Melia, F. 2006, ApJ, 636, L33 

Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 
Wachter, K., Leach, R., & Kellogg, E. 1979, ApJ, 230, 274 
Wang, Q. D., Lu, F. J., & Gotthelf, E. V. 2006, MNRAS, 367, 937 
Yusef-Zadeh, F, Bushouse, H., Dowell, C. D., et al. 2006, ApJ, 644, 198 
Yusef-Zadeh, F, Choate, D., & Cotton, W. 1999, ApJ, 518, L33 
Yusef-Zadeh, F, Melia, R, & Wardle, M. 2000, Science, 287, 85 
Zhao, J.-H., Young, K. H., Herrnstein, R. M., et al. 2003, ApJ, 586, L29 



References 



Anantharamaiah, K. R., Pedlar, A., & Goss, W. M. 1999, in Astronomical 
Society of the Pacific Conference Series, Vol. 186, The Central Parsecs of 



D. Porquet et al.: X-ray hiccups from Sgr A* observed by XMM-Newton, Online Material p 1 



Online Material 



D. Porquet et al.: X-ray hiccups from Sgr A* observed by XMM-Newton, Online Material p 2 



Appendix A: Absolute astrometry of the bright 
X-ray flare 

To obtain an absolute astrometry of the bright X-ray flare, we 
registered the positions of the soft X-ray sources in the third ob- 
servations on the Hipparcos celestial coordinate syste m using the 
X-ray counterparts of the Tycho-2 catalog's sources (H 0g et alj 
2000). We matched the positions of the pn sources detected 
in the 0.5-1.5 keV energy band obtained with the SAS task 
edetect_chain with the positions of Tycho-2 catalog's sources 
for the epoch of our observation. We found 4 X-ray counterparts 
of Tycho-2 sources using a cross-correlation radius of 5". The 
X-ray counterpart of Tyc6840 020 1, lying on a pn CCD gap, 
is not considered in the following. The offsets between pn and 
Tycho-2 positions of these three reference sources shows a sys- 
tematic translation. The weighted mean offset is -0V7 and 0'.'6 in 
right ascension and declination, respectively. Therefore, we cor- 
rect the X-ray position by subtracting this weighted mean offset. 
The residual RMS scatter in the corrected X-ray positions of the 
reference sources, U.'9, gives the astrometric uncertainty of the 
registration. This value is combined with the positional error of 
the PSF fitting to derived the final positional uncertainty. We list 
in Table lA.fl the corrected X-ray positions of the Tycho-2 refer- 
ence stars, the final positional uncertainties, and the distance to 
the Tycho-2 optical positions. The absolute position of the X-ray 

bright flare is aj 2 ooo = 17 h 45 m 40.0 s , 6 nm = -29°00'28.6" with 
a one-sigma positional uncertainty of l'.'O. 



Table A.l. X-ray counterparts of Tycho-2 reference stars. 



Name 


J2000 coordinates* 3 ' 


Err. (b) 


Name 


Dist. 


XMMU J 1745 


17 h 45 m 




Tyc 6840 


It 


25.7-285627 


25.7 s -28°56'27.8 


1.3 


666 1 


0.9 


43.8-291317 


43.8 s -29° 13' 17.7" 


1.7 


334 1 


0.9 


43.9-290456 


43.9 s -29°04'56.5" 


1.1 


590 1 


0.2 



(a) Absolute position corrected from systematic translation. 

<b) Final positional uncertainty obtained by combining the positional 
error of the PSF fitting with the astrometric uncertainty of the reg- 
istration on Tycho-2 sources. 



Appendix B: x 1 statistic versus W statistic 

The flare spectra are rebinned from 1 keV using grppha to have 
at least 25 counts per spectral bin. We ignore data above lOkeV, 
as for the fits with t he W statistic (see 93.21 ). XSPEC (version 
12.4.0;|Arnaud||l996) was used to fit the background-subtracted 
spectrum. For the October 2002, as for the W statistic spectral 
fit, we used as background the longer February 2002 observation 
(orbit 406, ObsID: 01 1 1350101) where no Sgr A* flare were ob- 
served. If the non-flaring spectrum observed on October 2002 
before the flare is used as background, very similar parameter 
fits are obtained. The errors and upper limits quoted correspond 
to 90% confidence level for one interesting parameter. The pa- 
rameter fits (see Table IB. IT ) are similar to that found using the W 
statistic for bright flares. For lower S/N flares the central values 
are different but consistent within the error bars with that the val- 
ues found using W statistic. The error bars are generally larger 
using the^ 2 statistic. 



Table B.l. Best fit parameters (using x 2 statistic) of the 
EPIC flare spectra for absorbed pegged power-law (pow), 
bremsstrahlung (brems), and black-body (bb) models, taking 
into account dust scattering (assuming A v =25mag). The errors 
are given at the 90% confidence level. 



Flare 


Model 




r/kT (b) 


r,roean(c) 
/Y 2-10keV 


^ 2 /d.o.f. <2 <d) 


April 4, 2007 


#2 


pow 
brems 
bb 


12.8+fj 
11 3 +1 - 9 

U,J -1,6 
7 q + 1.6 
' -1.3 


9 q+0.4 
z " J -0.4 

1-51! 


15.7+15 

i3.4 : : 

Q S+ - 8 
y - -0.7 


70.2/75 63.6% 
69.8/75 64.9% 
71.5/75 59.2% 


#3+#4+#5 (e) 


pow 
brems 
bb 


63«* 
6.4+- 

9 Q+2.6 


1 2+ - 7 

> 11 

9 i +0.6 


4 r 13 

A 1+ 6 - 6 
^• X -1.4 

3 3+ 12 

->">-0.7 


29.9/32 57.3% 
29.9/32 57.4% 
29.7/32 58.2% 


October 3, 2002 




pow 
brems 
bb 


12.113 
10.62! 
6.5+|;i 


9 9+0.3 
i.i„.3 
7 Q+3-4 

'• y -is 
1 6+° 1 


24 1+ 36 
20.9+^ 


90.4/102 78.8% 
90.9/102 77.6% 
98.1/102 59.0% 



<il) N H is units of 10 22 cm 2 . 

<b) The black-body temperature is given in keV. 

<c) Mean unabsorbed fluxes for the flare period in the 2-10 keV energy 
range in units of 10~ 12 erg cm 2 s _I . 

<d) Q is the null hypothesis probability, i.e., the probability of getting a 
value of x 1 as large or larger than observed if the model is correct. 

(e) These parameter fits for flare #3+#4+#5 are only given to illus- 
trate the impact on both the central parameter fit values and the 
constraints on error bars, when using x 1 statistics in case of low 
statistics, instead of W statistics (see Table[3]for comparison). 



Appendix C: Reprocessed light curves and spectra 
of the X-ray flares previously observed with 
XMM-Newton 

We report the reprocessed light curves (Fig. IC.U and spectra 
(Fig. IC.2l and lC.3b of the X-ray flares observed on October 2002 
and M arch 2004 that were report ed previously by Porquet ~et al.l 
(120031) and lBelanger etalJ (120051) . respectively. Table ICTl gives 
for each flare: the start and end times, the duration, the to- 
tal EPIC quiescent-subtracted counts, and EPIC quiescent- 
subtracted count rate at the peak flare. 

The X-ray flare observed on March 2004 was contaminated 
by a high level of background proton flares which lead for pn 
to numerous switches between science and counting mode, and 
therefore to a significant loss of exposure time. Consequently, a 
bin time interval of 100 s contains less than 50% of (scientific 
mode) exposure, which is not enough to obtain a good estimate 
of the count rate. We used a bin time interval of 150 s, which is 
sufficient to obtain a good estimate of the count rate, when com- 
paring with the MOS light curves. Obviously, this loss of expo- 
sure on this moderate flare affected the quality of its pn spectrum, 
and the resulting constraints derived on its physical parameters. 
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Table C.l. Characteristics of the X-ray flares observed in October 2002 and March 2004 (see Fig. IC.lt . 



Flare 


Start time <a) 
hh:mm:ss s 


End time (a) 
hh:mm:ss s 


Duration 

s 


Total (b) 
cts 


Peak (c) 
cts s" 1 


Det. (d) 

(T 


. peak (e) 
^2-10 keV 

10 34 erg s _1 


Ampl. (f) 


03/10/2002 
31/03/2004 


10:07:40 
22:52:24 


150026860 
197160744 


10:54:20 
23:12:24 


150029660 
197161944 


2800 
1200 


2158 
447 


1.510 
0.653 


29.6 
12.8 


35.9^ 
11.0™ 


(g) 149 +21 

46 +9 



Start and end times of the flare time interval defined as the period where the EPIC light curve deviates from the non-flaring level at a confidence 
of 95%. See CTilfor details. 

Total EPIC counts in the 2-10 keV energy band obtained during the flare interval after subtraction of the non-flaring level. 
EPIC count rate in the 2-10 keV energy band at the flare peak after subtraction of the non-flaring level. 
Detection level at the flare peak in a. 

2-10 keV luminosity at the flare peak, assuming an absorbed power law model (taking account dust scatter effect), see parameter fits in Table[3] 
Amplitude of the flare defined as the ratio of the 2-10 keV flare lu minosity at the peak and the 2-10 keV quiescent luminosity of Sgr A* 

o bserved with Chandra (i.e., 2.4 x 10 33 erg s -1 , Baganoff et al. 2003, and this work, see <j4j. 

In Porquet et al. (2003), the amplitude was calculated using the 2-10 keV quiescent luminosity of Sgr A* reported in Baganoff e t al.l ( 120011) . 
i.e. 2.2 x 10 33 erg s" 1 . Using this latter quiescent luminosity value, we would find for October 2002, an amplitude of 163+j 3 . 
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Fig. C.l. XMM-Newton EPIC spectra of the flare observed on 
October 2002. The spectrum for the non-flaring period with the 
same extraction region is used as background. Squares, triangles 
and stars indicate MOS1 (red), MOS2 (green), and pn (black) 
data, respectively. The lines show the best-fit model using an 
absorbed power-law model taking into account dust scattering 
(see Table IB. 11 for best fit parameter values; pow model). The fit 
residuals are shown in the bottom panels. 
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Fig. C.l. Light curves of reprocessed X-ray anterior flares ob- ™x 
served with XMM-Newton. Left: October 2002 flare. Right: < 
March 2004 flare. -2 
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Fig. C.3. Same as Fig. lC.2l but for the flare observed on March 
2004. 
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